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Whereas it is easy to reduce the translational symmetry of a molecular system by 
using, e.g., Jacobi coordinates the situation is much more involved for the rotational 
symmetry. In this paper we address the latter problem using holonomy reduction. We 
suggest that the configuration space may be considered as the reduced holonomy bundle 
with a connection induced by the mechanical connection. Using the fact that for the 
special case of the three-body problem, the holonomy group is 50(2) (as opposed to 
50(3) like in systems with more than three bodies) we obtain a holonomy reduced con- 
figuration space of topology R+ x 5 1 . The dynamics then takes place on the cotangent 
bundle over the holonomy reduced configuration space. On this phase space there is 
an 5 1 symmetry action coming from the conserved reduced angular momentum which 
can be reduced using the standard symplectic reduction method. Using a theorem by 
Arnold it follows that the resulting symmetry reduced phase space is again a natural 
mechanical phase space, i.e. a cotangent bundle. This is different from what is ob- 
tained from the usual approach where symplectic reduction is used from the outset. 
This difference is discussed in some detail, and a connection between the reduced dy- 
namics of a triatomic molecule and the motion of a charged particle in a magnetic field 
is established. 
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1 Introduction 



In molecular dynamics which is the subject of this paper and generally in dynamical sys- 
tems theory the reduction of the number of degrees of freedoms is of central importance for 
both computational and conceptual reasons. A molecular system is a many body system 
consisting of the nuclei and electrons of the constituting atoms. The electronic degrees of 
freedoms are typically dealt with in a Born-Oppenheimer approximation. Since the nuclear 
masses are a few thousands times bigger than the mass of an electron one assumes that the 
nuclei adiabatically interact via the forces obtained from a potential energy surface that is 
obtained from the electronic ground state energy as a function of the nuclear configurations. 
The computation of such potential energy surfaces is based on density functional theory and 
other methods and is an art in physical chemistry. For several molecular systems such poten- 
tial energy surfaces are tabulated in the chemistry literature. Given such a potential energy 
surface a molecular system reduces to an iV-body system which only involves the degrees of 
freedom of the N nuclei in the system. This iV-body system can then be treated classically 
or quantum mechanically. In particular for light atoms (respectively nuclei) like hydrogen 
quantum effects might play an important role which make a quantum mechanical treatment 
necessary. We note that state of the art quantum computations for, e.g., reactive scattering 
are even today only feasible for three or maximally four atoms. For this reason and also 
conceptual reasons one desires to get rid of as many degrees as possible. A reduction of the 
number of 'effective' degrees of freedom of a molecular iV-body system can be achieved by 
exploiting the symmetries of the system. These symmetries consist of overall translations and 
rotations. The reduction of translational degrees of freedom is simple and can be achieved 
by using Jacobi coordinates or changing to a centre of mass coordinate system. For rota- 
tions, the situation is much more involved as a clear distinction between rotational degrees of 
freedom and (internal) vibrational degrees of freedom only exists in an approximate sense in 
the vicinity of an equilibrium position. Here the distinction between vibrations and rotations 
can be achieved from the so called Eckart frame ]X\ that is widely used in applications [2] . 
This approximation is however only of local validity since large amplitude vibrations may 
produce rotations. A major step towards a geometric understanding of why a separation of 
rotations and vibrations cannot be achieved globally goes back to the work of Guichardet 
[3] who used the differential geometry framework of principal bundles to give a mathemat- 
ically rigorous definition of vibrational motions. He showed that the translational reduced 
configuration space is a principal bundle with structure group given by the special orthogo- 
nal group, and introduced a connection which naturally relates to molecular motions. The 
inseparability of rotations and vibrations then follows from the nonvanishing curvature of 
this so called mechanical connection. Iwai and Tachibana [U E] used Guichardet 's approach 
to study in great detail both the classical and the quantum mechanical dynamics of iV-body 
molecular systems. Using the setting of principle bundles Iwai [6] in particular showed that 
the Eckart frame can also be defined for general configurations (i.e., no necessarily equilib- 
rium configurations) of a molecule. However, this frame is then not unique and therefore 
not suitable for studying large amplitude vibrational motions of a molecule. Iwai moreover 
applied the Marsden-Weinstein-Meyer symplectic reduction procedure [TJ [S] to reduce the 
constant angular momentum motion of an iV-body system. He showed that for nonvan- 
ishing angular momentum the reduced phase space is then no longer a natural mechanical 
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system in the sense that it is no longer given as the cotangent bundle over a (reduced) con- 
figuration space. A gauge theoretical interpretation of the reduction of symmetries and the 
related choice of a reference frame in iV-body systems was introduced in [21 [TO]- In their 
constructive and instructive paper Littlejohn and Reinsch[10j used Lagrangian reduction in- 
stead of symplectic reduction mentioned above. For more related work we mention the refer 



In this paper we use modern tools from the geometric description of molecular motion 
described above to introduce a new way to reduce the symmetry specifically of triatomic 
molecular system. We obtain a reduced configuration space and deduce the reduced dynamics 
for a triatomic molecule in a way which can be summarized as follows. Consider three atoms 
(or nuclei) in R 3 . Using Jacobi coordinates the translational symmetry in the absence of 
external forces can be used to reduce the nine- dimensional configuration space R 3 x R 3 x R 3 
of the triatomic system to the six-dimensional space R 3 x R 3 . Excluding collinear (and hence 
also collisional) configurations from R 3 x R 3 gives the translation reduced configuration space 
P on which the special orthogonal group 5*0(3) acts freely. The space P is a principal bundle 
with base space given by the positive half space R^ [TTJ. Kinetic energy gives a metric on P, 
and a connection can be obtained by defining horizontal spaces as orthogonal complements 
of the tangent spaces of orbits of the 5*0(3) action. As known [3J [TD] the connection on P has 
a nontrivial holonomy group which is 50(2). This enables us to use the holonomy reduction 
theorem [19] to reduce P to the holonomy bundle which we denote by Q. Since P is a trivial 
bundle [TTJ, Q is also trivial and hence topologically R+ x 50(2), or equivalently R 3 x 5 1 . 
The reduced phase space is then given by the cotangent bundle T*Q. We explicitly derive 
the Hamiltonian on T*Q and deduce the reduced dynamics on T*Q. In the final step we then 
use the conservation of the reduced angular momentum related to an 5 1 action on T*Q to 
apply the symplectic reduction procedure. Using a theorem in [17] we find that the reduced 
phase space is then a natural mechanical system, namely the cotangent bundle over Q/S 1 . 

We note that there is no natural way to generalize these results to systems of four or more 
atoms. The reason is that triatomic systems are in many respects special. For example, the 
holonomy group of a system of four or more atoms is 50(3), and the translation reduced 
space is not a trivial bundle [TO] . 

2 Reduced configuration space 
2.1 Principle bundle picture 

Consider a molecular system of three atoms. Let Xj G R, % = 1, 2, 3, be the position vectors 
of these atoms. Suppose that there are no external forces. Then the mass-weighted Jacobi 
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can be chosen to reduce the symmetry of overall translations. (For different choices of Jacobi 
vectors see[Bj) Excluding collinear (and hence also collisional) configurations we obtain the 
six-dimensional translation reduced configuration space 

P = [ x = ( r ,s) : Ar + /is ^ for all (A,/i) G R 2 \{0}} C R 3 xR 3 . 

Proper rotations g G 50(3) act on P in the natural way 

#(r,s) = (gr,gs) . 

On P this action is free and it thus follows from standard results that 

M := P/S0(3) 

has a manifold structure. The space M is usually referred to as shape space or internal space. 
Furthermore, the canonical projection 71 : P — > M defines a principal bundle with structure 
group 50(3) [3J. This means that P consists of smoothly glued copies of 50(3), i.e., locally, 
P is diffeomorphic to M x 50(3). Topologically, this local decomposition also holds globally 
which following Iwai [11] can be seen as follows. Using Jacobi coordinates 



r = a/ (r,r), s = a/ (s,s), = cos 1 ((r, s) /rs) , 

where (•, •) is the usual dot product on R 3 , and introducing coordinates 

W\ = r 2 — s 2 , w 2 = 2rscos(j), w 3 = 2rssin0 > 

one sees that M = R 3 = {(wi, w 2 , w 3 ) : w 3 > 0}. As pointed out in [11], P is a trivial bundle 
as M is contractible to a sine le point. So, topologically, P = R 3 x 50(3). 

2.2 Nontrivial holonomy 

Turning back to the action of 50(3) on P one can see that the fundamental vector field A 
associated with an element A in the Lie algebra so(3) is given by 



~ d 
Alx= Jt 



(e tA x) , (1) 

t=o 



or equivalently, 

A\ x = (At, As) = (w x r, w x s), (2) 
where w G R 3 is the unique vector corresponding to A by the isomorphism 

/ -a 3 a 2 \ I at \ 
R- 1 : so(3) -> R 3 , a 3 -ai a 2 . (3) 
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Let be an orbit of the 50(3) action, say iV = 50(3)x for a point x G P, then T x .iV = 
: A G so(3)|. Consider the orthogonal complement H x of T X N in T X .P with respect to 
the Euclidean dot product on P given by 

dx 2 = (r, r) + (s, s) . (4) 
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Clearly the distribution x t— > H x , which we call the horizontal distribution, defines a connec- 
tion [IT] uj : TP — > so(3) on P which is a special case of the mechanical connection defined 
in [IS]. A vector field X* with X*\ x e H x for all x G P is called horizontal. The horizontal 
lift of a vector field X on M is accordingly the unique horizontal vector field X* on P such 
that dir(X*) = X. We have u(X*) = for every horizontal vector field X* and 00(A) = A for 
every fundamental vector field A. In order to compute the horizontal lifts of the coordinate 
vector fields d r ,d s , on M we give an explicit expression for the metric dx 2 in 01]). To this 
end we follow [TQJ [T3] and introduce a frame u l5 u 2 , u 3 in R 3 according to 

r = rui, 

s = s cos^U! + ssin0u 2 , 
u 3 = Ui x u 2 . 

If Euler angles (a, (3, 7) on £0(3) are chosen via 

g = e ^(aei) e ^e 2 ) eJ R( 7 e 3 ) ) < «, 7 < 2vT, < /3 < 7T, 

where ei, e 2 , e 3 is the standard basis of R 3 , R is defined in ([3]) and g>ej = Uj, i = 1, 2, 3, then 
with 

©i = sin. j d(3 — sm(3 cos ^ da, 

9 2 = cos •yd/3 + sin (3 sin 7 da, 

9 3 = cos j3 da + d'j, 



one obtains [13 



dr = dr u x + r0 3 u 2 — r0 2 u 3 , 

(is = ?7iUi + 7/ 2 U 2 + T] 3 U 3 , 



where 



r\i = ds cos <f) — ssinipd(f) — ssin00 3 , 
7] 2 = ds sin cj) + s cos 4>d(f> + s cos 4> 6 3 , 
r/ 3 = ssin^Oi — scos</>0 2 . 

In local coordinates the metric <ix 2 then assumes the form 

dx 2 = dr 2 + r 2 (6 2 + 6 2 ) + r£ + rf 2 + V 2 3 . 

This expression can be used to locally compute the horizontal lift X* of a vector field X on 
M: X* is orthogonal to 8^8^,8^, and dn(X*) = X. It follows that 

8* r = 8 r , d* s = d s , dl = d^- J 2 <9 7 . (5) 

T ~\~ S 

In gauge theory the factor 2 r ^_ 2 is referred to as a component of a Yang-Mills potential [5] . 
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2.3 Holonomy reduction 

By Equation ([5]) we have arrived at the well-known phenomena of inseparability of rotations 
and vibrations |3j [6]. Namely from (J5j we see that the distribution spanned by d*,d*,d^ 
is not integrable, and hence, if these vector fields are considered as infinitesimal vibrational 
motions one can say that vibrations generate rotations. This is why the internal space M is 
not a submanifold of P J5J. On the other hand d*,d*, d% d 1 do span an involutive and hence 
integrable distribution. The maximal integral manifold Q x of that distribution at a point 
x G P is a good candidate for being the reduced configuration space because vibrational 
motions through x live in that space. In fact we will obtain the reduced dynamics of a 
triatomic molecule on the cotangent bundle over Q x by employing the holonomy reduction 
of principle bundles: A curve on P is called horizontal if its tangents are horizontal. Fix 
a point x G P and denote by P(x) the set of all points in P which can be joined to x by 
horizontal curves. It is known that [31 [10] the holonomy group of u is 50(2) (see also [A]) . 
and since M is connected and paracompact the holonomy reduction theorem [19] implies that 
P(x) is a reduced bundle with structure group 50(2), which is in fact Q x . Furthermore, Q x 
is a trivial bundle as it has the same base space as P. These observations suggest that the 
reduced configuration space of a triatomic molecular system is topologically x 50(2). 
The induced metric on Q x is thus 

^ = * 2 + * 2 + ^^ 2 + ^C 2 , (6) 

where 

( = s 2 d(j) + (r 2 + s 2 )d-f. (7) 

3 Reduced dynamics 
3.1 Angular momentum 

In the following we want to put our derivation above into the context of some well known 
results. It is known [1] that in the case of small vibrations one can separate vibrations and 
rotations in the vicinity of an equilibrium point. In the present situation if one chooses 
g?7 = in (jEJ) the well-known Eckart kinetic energy is obtained. This is the gauge dependent 
internal metric h^ u in [10]. Thus one can conclude that in case of small vibrations the internal 
motions of molecule live in the integral manifolds of the distribution spanned by d r , d s , d^, 
called the Eckart space. Next, consider the angular momentum 

J = R(r x dr + s x ds) 

on P which is computed locally to be 

J = R((r 2 Q 2 — s cos0?73)u 2 + (r 2 Q 3 + s cos0r] 2 — ssin07/i + s sin (f>r] 3 )u 3 ) . 

So, its restriction to Q x is 

3\ Qx = R((u 3 ). 
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If the angular momentum of the system is identically zero, then ( = 0, and hence 



2„2 



dq 2 = dr 2 + ds 2 + ' 



r z + s z 

3.2 Holonomy reduced Hamiltonian 

In the case of vanishing angular momentum the Hamiltonian is obtained to be 

H = -p 2 r + \pl + \{\ + \)pl + V(r, s, 0), 



2 V 2 



where p r ,Ps,P<j> are the conjugate momenta and V(r, s,0) is the potential energy which is 
assumed to be rotationally invariant. This Hamiltonian is widely used in applications. By 
(J7J) we observe that vibrational motions live in the integral manifold of the distribution 
spanned by <9*,<9*,<9 7 . That space may be called zero-angular momentum space. 

If the reduced angular momentum J\q x is a non-zero constant, we have ( = const. ^ 0. 
Then, equivalently, is a non-zero constant and hence the vibrational motions remain in a 
three-dimensional affme space which is parallel to the zero-angular momentum space. 

Taking into account the contribution of ( in the induced metric dq on Q x in ([6]) the 
Hamiltonian in ([8]) changes to the general holonomy reduced Hamiltonian 

11111 1 1 

H = -P 2 r + -P 2 S + -(^ + -^)pl - ^P^ + + V{r, s, 0). 

The corresponding Hamiltonian vector field is given by 

X = p r d T + p s d s + + -^)^ - ^P-}j d<t> + -^(p-y - P4,)d 7 + 

Since 7 is cyclic the conjugate momentum p 1 is conserved. To put it another way J| Q ^ is 
an S' 1 -equivariant momentum and the standard symplectic reduction theorem can be applied. 
Using a theorem by Arnold (see [T7], page 378) the resulting reduced phase space is again a 
natural mechanical system, i.e. a cotangent bundle. 



4 Comments on related work 

4.1 The relation between the motions of a triatomic molecule and 
a charged particle in a magnetic field 

In [201 EI] the idea is introduced to describe the motion of a charged particle in a magnetic 
field by extending the configuration space R 3 to R 3 x S 1 such that the angle corresponding 
to S 1 is cyclic and its conserved conjugate momentum gives the charge of the particle in 
the magnetic field. Since the holonomy reduced configuration space is R 3 ^ x S* 1 we can 
identify the motion of a triatomic molecule to that of a charged particle in a magnetic field 
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as follows. Let q denote a point in M with coordinates (r, s,0). If A denote the one-form 
(r 2 + s 2 ) _1 C = (r2 ! ) _ s 2) ^0 + d'j on R/j_, then by the metric the kinetic energy can be written 
as 

Lk = \\H? + ^(A-q + 7) 2 

which is reminiscent of the so called Kaluza-Klein Lagrangian [21] • The conjugate momenta 
are then 

P = -^r- = q + (A-q + 7 )A 

and 

dL K 

Vi = ~q7- = A • q + 7. 

The one-form A plays the role of a vector potential for the magnetic field. The conserved 
momentum p 1 is the charge e = cp 7 (with c denoting the speed of light) [2T| . 



4.2 Relation to symplectic and dimensional reduction 

In [6] the symplectic reduction procedure [3, [E] is applied to the N— body problem. The 
cotangent bundle over the translation reduced configuration space P is a symplectic manifold 
with the canonical two-form, and the angular momentum J : T*P — > so(3) is an equivariant 
momentum map. For a fi 7^ 0, it is shown that J~ 1 (fi) is a principal bundle with structure 
group SO (2) whereas the zero momentum space J _1 (0) is a principal bundle with structure 
group SO (3). Furthermore J~ 1 (0)/S'O(3) is shown to be diffeomorphic to T*(P/SO(3)), but 
J~ l (lj) / SO{2) is no more a cotangent bundle because of dimensionality. As pointed out in 
[6] the procedure for the latter when applied to for three-bodies is in fact the elimination of 
nodes. 

In contrast to the symplectic reduction procedure the first step in this paper was to pass 
from the translation reduced configuration space P of a triatomic molecule to a subbundle 
Q (the holonomy reduced bundle) which is a principle bundle with structure group 5*0(2). 
Afterwards the angular momentum is then restricted to T*Q, and finally the sypmlectic 
reduction procedure is applied. The reduced space is then always a cotangent bundle as 
follows from a theorem by Arnold (see [T7], page 378). 

We note that the method used in the present work is strongly related to dimensional 
reduction [22} [23], a method developed for symmetries of gauge fields. More precisely, in the 
case of spherical symmetry in 6 dimensions applied to an SU(3) gauge theory, the 2 extra 
dimensions describing a sphere of radius R. One solution, with the largest set of Higgs fields, 
reduces to the 4-dimensional Weinberg-Salam model without fermions [22]. 



5 Conclusions 

In this paper we used the geometric theory of molecular mechanics [3j [HI EE] to reduce the 
number of degrees of freedom in the molecular three-body problem. We followed the principal 
bundle setting of Guichardet [3] on the translation reduced configuration space, and using the 
holonomy reduction theorem [TH] it was possible to reduce to a principle subbundle. This 
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may be interpreted as separating two rotational degrees of freedom from the maximal space 
that includes vibrational motions. It was then possible to induce the angular momentum and 
apply the very symplectic reduction procedure (to be precise, we used it in the form of the 
Noether's theorem here). This way, the remaining momentum space which is of 6 dimensions 
and also a phase space was obtained. The computations were local for the purpose illustration 
but the method is intrinsic. In some sense, the resulting space which is of 3 dimensions can be 
seen as a rotationless space. For the case of zero angular momentum this is a known earlier 
result. In the present paper it was generalized to the case of non-zero angular momentum. 
In particular we used our approach to rephrase the well known fact [20] that a triatomic 
molecular system behaves as a single particle in a magnetic field. 

A A lemma by Guichardet 

For completeness, we give a brief proof of the fact in that a vibrational motion of a triatomic 
molecule, which is defined as a curve with horizontal tangents, remains in a fixed plane as 
originally formulated by Guichardet [3J: Let x(t) = (r(i),s(i)) be a horizontal curve on P. 
We show that F x ^ : = span {r(£), s(£)} is fixed. Indeed, since x(t) is horizontal it is orthogonal 
to all fundamental vector fields which are given in (|2J), and hence r(t) x r(t) +s(t) x s(t) = 0. 
Let y(t) be a curve in R 3 with (y(t),y(t)) = 1 which is orthogonal to F x ( t y So, y(i) is 
orthogonal to F x{t) . Hence (y(t),r(t)) = (y(t),s(t)) = so (y(t),r(t)) = (y(*),s(t)) = 
which implies y is zero. 

As a conclusion of the above fact it is observed [10] that during vibrational motions or 
shape deformations the Jacobi vectors remain in a fixed plane, and hence the Jacobi vectors 
before and after the vibrational motion can be transformed to one another by a plane rotation, 
i.e. the holonomy group is £0(2). 

B The kinematic group 

Different clusterings of position vectors give rise to different choices of mass-weighted Jacobi 
vectors. These different choices are related to each other by transformations which are called 
democracy transformations [10J. The set of all democracy transformations forms a subgroup 
of the symmetry group 5*0(3) called democracy or kinematic group. For the three-body 
problem the kinematic group is 5*0(2). This is another special feature of the three-body 
problem. 
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